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Abstract 

The recent measurement from the CLEO experiment presents the DD*n coupling, 17.9±0.3±1.9. 
This value is much larger than any of QCD sum rule predictions available in literature. We report 

o, 

that, with a relevant treatment of the continuum subtraction as well as with the asymptotic form of 
' the twist- 2 pion wave function, the light-cone QCD sum rule can provide the coupling comparable 

m ' 



to the experimental value. The stability of the resulting sum rule becomes much better with these 
corrections. 
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I. INTRODUCTION 



Recently, the CLEO collaboration presents a new measurement of the D*Dir coupling, 
9d* Dn — 17.9 ± 0.3 ± 1.9. Though this value is somewhat consistent with quark model 
predictions [0 and lattice calculations 0, it is much larger than any of QCD sum rule 
predictions, which gives some skepticism in accepting the various QCD sum rule results. In 
particular, the conventional QCD sum rules relying on the short- distance expansions |f]. [5|. 

[7J consistently restrict the coupling go*Dn < 10. In these sum rules, QCD inputs are 
well-constrained by the low-energy theorem. Thus, the current disagreement may due to 
the fact that various pion matrix elements appearing in the conventional sum rules may 
not converge fast enough in the short-distance expansions. One may need to do a partial 
summation to all order of those matrix elements. 

The other method is to construct QCD sum rules from three-point function ||, which 
leads to the coupling around 6, again much smaller than the experimental value. Its recent 
extension ||] using exponential type parametrization for the form factor seems to give a 
much larger value around 15. However, it should be remembered that the operator product 
expansion (OPE) of a three-point function leads to an unphysical behavior at high momen- 
tum transfers JT^|, which in fact was one of the motivation for constructing light-cone QCD 
sum rules (LCQSR) relying on the expansions along the light-cone. Also, LCQSR may im- 
prove the conventional QCD sum rules by having the pion wave functions that encode the 
partial summation to all order of the pion matrix elements. However, LCQSR [O, [12] yields 
the D*Dtt coupling only about 12, still smaller than the experimental value. Thus, as far as 
the D*Dix coupling is concerned, all the sum rules do not provide the experimental value. 
Furthermore as quark model and lattice calculations [§ agree with the experiment, it is 
important to re-investigate the existing sum rules and look for a way to reach an agreement 
with them. 

One crucial input in the prediction of LCQSR is the twist-2 pion wave function at the mid- 
dle point ^^(0.5) . Refs.0, Q in their sum rules use ^(0.5) ~ 1.2 obtained by comparing 



two different Lorentz structures of the light-cone sum rule for the pion-nucleon coupling [|T3" 



However, as the pion-nucleon sum rule has a strong dependence on the Lorentz structure |14j 



considered, ^,,-(0.5) ~ 1.2 is questionable. Thus, one may attribute the current discrepancy 
in the D*Dn coupling to the uncertainty of (p n (0.5). Indeed, there are some suggestions that 
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the twist-2 pion wave function should be of the asymptotic form whose value at the middle 



point is ^(0.5) ~ 1.5 JT5 , [T6[1 . As we will discuss below, this new input only increases the 
D*Dn coupling by 9 %, still not enough to reproduce the experimental value. Certainly, the 
existing sum rules need further improvements. 

In this letter, we suggest one possibility to improve the existing LCQSR. Specifically, we 
show that the continuum subtraction in the existing LCQSR is mathematically ill-defined. 
A similar suggestion was reported in Refs. [[[]], [18], In particular, we demonstrate 

that the OPE of the form J Q du/[m 2 c — up\ — (1 — w)pf] should not entail the continuum 
subtraction if one strictly follows the QCD duality assumption in constructing the continuum 
contributions. This correction combined with the asymptotic twist-2 pion wave function 
leads to the D*D7r coupling comparable to its experimental value. The resulting sum rule 
is stable against the variation of the Borel mass. 

II. THE CONTINUUM CONSTRUCTION 

In this section, we demonstrate a possible modification in obtaining the continuum sub- 
traction in the existing LCQSR [|ll|, [L2| . The leading OPE in the existing LCQSR takes the 
form 

dU (1) 



m 2 — up{ — ( 1 — u)p2 

where m c the charm quark mass and P\,P2 are the two momenta associated with the coupling 



vertex. The leading OPE in Refs. [JTT], |12| also contains the pion wave function in the 
numerator of the integral but here we suppress that for a mathematical simplicity. The 
main issue regarding the continuum subtraction is not affected by this simplicity. 

It was claimed in Ref. |TT| |12|| that, after the continuum subtraction and the double Borel 
transformation, the OPE, Eq.(|l|), appears in the final sum rule as the factor 

M 2( e -ml/M* _ e -S /M^ (2) 

where M 2 is the reduced mass of the two Borel masses associated with the double Borel 
transformations. We will show that the continuum subtraction factor, M 2 e~ So ^ M2 , comes 
from a mathematically spurious term and it should not be a part of the final sum rule. 
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To determine the continuum subtraction for the given OPE, one needs to determine the 
spectral density p(si,s 2 ) firstly from the double dispersion relation, 



L 



-r iSi r iS2 i eM _ . (3) 



, m? c - up\ - (1 - u)p 2 2 J J (si - pf)(s 2 - vl) 
Note that the lower boundaries of the integrals in the right-hand side should be chosen such 
a way that the entire region of p(s±, s 2 ) 7^ is included. According to the prescription given 



in Ref . [fUJI , the spectral density is given by 

p{s 1 ,s 2 ) = 6(sx - s 2 )9(s 1 - m 2 c ). (4) 

An alternative way to obtain this spectral density is to take the imaginary part of Eq.(|T]) 
with respect to the two momenta, p\ and p\. 

A common practice in LCQSR is to (1) put this spectral density back to the double 
dispersion integral, (2) restrict the integral below the continuum threshold So (according 
to QCD duality), J °° ds\ J °° ds 2 — > J Q So dsi f Q So ds 2 , and (3) take the double Borel trans- 
formations. This procedure as advocated by Belyaev et.al. precisely yields the simple 
prescription given in Eq. ([|). We want to point out that the second step is dangerous as it 
yields a spurious contribution. To show this, let us see how the spectral density Eq.(§) re- 
produces the OPE Eq.(jl]) within the double dispersion relation Eq.(|3|). We put the spectral 
density Eq.(f|) in the double dispersion relation and obtain 

dsi I ds 2 - ^ 2T = f du [ dsi- — c ^ 2 . (5) 



^0 



si -Pi)(s 2 -Pi) Jo Jo (si - upj - (1 - u)p 2 2 ) 2 



Note, we have used the Feynman parameterization to get an integral over the Feynman 
parameter u. We then perform the integration by part to separate into the two terms 



1 '- x S( Sl -m 2 c ) f 1 . 0(si-m 2 c ) 



du / dsi 2 77^ — r~2 — / du- , 

si - up{ - (1 - u)pi Jo Si - up( - (1 - u)pt 



si=oo 



(6) 

si=0 



10 JO °1 ~ u>V\ ~ V s - ~ lL )f2 JO °1 "Fl — V x — u Jf2 

One can see that the first term precisely yields the anticipated OPE while the second term, 
though vanishes under the present boundaries, is an unnecessary spurious term. When 
QCD duality is applied (i.e. restricting the integral below the continuum threshold So), it is 
mathematically sensible to apply to the first term only. Any contribution from the second 
term is ill-defined as the second term is mathematically spurious. However one can show that 
the continuum subtraction factor M 2 e _5o//A/2 comes from the spurious second term. When 
the integral interval in Eq. (||) is switched to ~ 5*0, the second term becomes du/[So — 



upl — (1 — u)p2], which no longer vanishes. Then under the double Borel transformations 
the second term precisely yields the continuum subtraction factor M 2 e~ 5o//M2 . Therefore, 
the continuum subtraction factor M 2 e~ s °^ M2 should be dropped in the final expression of 
the sum rule as it is from the mathematically spurious term. 

Once the factor M 2 e~ s °^ M2 is dropped, then the resulting sum rule does not depend on the 
continuum threshold, indicating that there is no continuum contribution. Intuitively, absence 
of the continuum contribution may seem strange because the current in the correlation 
function can couple to higher resonances as well as the lowest resonance. But one can show 
that a s corrections to the perturbative part can give a continuum contribution within our 



prescription as the corrections are logarithmic functions of Si and S2 [ 20 1 . Therefore one can 
pick up small continuum contribution [at an order and, in the present calculation 

without a s correction, having no continuum contribution is not against the intuitive picture. 

III. RE-CALCULATION OF THE D*Dtt COUPLING 

Having suggested a modification in the previous LCQSR, we now re-analyse the sum rule 



for the D*Dn coupling. We take the sum rule formula for the D*Dtc coupling from Ref.flT 



and obtain the solid curve in Fig. [T]. This is the same curve that was presented in Ref.[]TT 



Based on this curve, Ref. JTTJ obtained JdJd*9d*d-k = 0.51 GeV 2 , which yields go'Dw = 12.5 



if we use fo = 170 MeV and fo* = 240 MeV obtained from two-point vacuum sum rules. 
The dashed curve is obtained when we take the asymptotic twist-2 pion wave function, 
ip n (0.5) = 1.5, while keeping all other parameters fixed as in Ref. ||1 1|| . This correction 
increases the coupling only by 9 %. Thus, with this correction only, LCQSR cannot still 
reproduce the experimental value of the coupling, 17.9. 

However the correction coming from the continuum factor is substantial. The dot-dashed 
curve is obtained when the continuum subtraction factor M 2 e~ s °/ M2 is taken out in the sum 
rule according to our suggestion in Sec.|I|. One clearly sees that this curve can provide a 
much larger coupling. Furthermore, the curve has a minimum around which the variation 
with respect to the Borel mass is minimized. This indicates that, around the minimum, the 
sum rule is well saturated by the OPE terms included in the calculation. For a qualitative 
estimate, instead of doing a detail analysis, we take the minimum of the curve to calculate the 
coupling. The minimum gives f d fD*9D* d-k = 0.66 GeV 2 , which in turn yields the coupling 
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9d*d-k = 16.2. This value is only 10 % smaller than the experimental one. The other source 



of error is the uncertainty in the decay constants fjj and fp*. Indeed, from Ref . ||1 If , the 
decay constants are known to be fu = 170 ± 10 MeV and fp* = 240 ± 20 MeV. If we take 
these uncertainties into account we obtain the range 

14.1 < g D *Dn < 18.75 (7) 

which certainly overlaps with the experimental value. 

Of course, there are other possibilities to improve the existing LCQSR. For example, the 
current LCQSR is up to twist-4. There might be some contributions from higher twists. Also 
the twist-4 element denoted by 5 2 is currently known to be 0.2 GeV 2 based on QCD sum rule 
calculation of Novikov et. al. [^1] This value however crucially depends on the factorization 



assumption for the four-quark condensate ((qq) 2 ). Based on the rho meson sum rule [2^], the 
four-quark condensate can be much larger, which may shift S 2 to a larger value. Nevertheless, 
what we want to emphasize is that all these improvements are expected to be less than 10 
% because, otherwise, the twist expansion can not be valid. To be comparable with the 
CLEO experiment, one needs a correction of order 40 % in the existing LCQSR calculation. 
Therefore, the continuum correction that we are addressing in this work provides the most 
important modification to the existing LCQSR and it helps to reach an agreement with the 
CLEO experiment. More interesting is that, with this improvement, LCQSR provides the 
coupling comparable to the quark model and lattice calculations M. 
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FIG. 1: fjj f d* 9 D* Dir versus the Borel mass in the light-cone QCD sum rule. The solid line is from 
the sum rule of Ref. [pTlJ] - The dash curve is obtained when (p w (0.5) = 1.5 is used. The dot-dashed 
curve is obtained with the continuum correction. 




S 



